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Motivated by the observation of localized circular excita- 
tions ('oscillons') in vertically vibrated granular layers Q, 
we numerically investigate an extension of a Swift-Hohenberg 
model that exhibits a subcritical transition to square pat- 
terns. For sufficient subcriticality, stable oscillon structures 
are found. The localization mechanism is quite general and 
is due to non-adiabatic effects. Bound structures of oscillons 
of equal and opposite polarity are found with bound states 
of like polarity being less robust. Much of the phenomena 
are consistent with the experimental observations and suggest 
that oscillons are not specific to patterns in granular media 
or to parametrically driven systems. Experimental tests are 
suggested that would determine whether this minimal frame- 
work is sufficient to describe the phenomena. 

PACS: 47.20.-k; 47.20.Ky; 47.54.+r; 81.05.Rm 
Keywords: Fronts; Granular media; Localized states; Para- 
metric excitation; Pattern formation; Square patterns 



I. INTRODUCTION 

Recent experiments on vertically vibrated thin layers 
of granular material have revealed a wealth of interesting 
phenomena |l],^,|],[|]||| ■ Depending on the driving fre- 
quency and amplitude different types of standing-wave 
patterns are observed. They can be ordered, e.g. stripe 
patterns and square patterns as they are also observed 
in parametrically excited surface waves on liquids, or 
disordered. Particularly exciting is the observation of 
two-dimensional localized structures which were termed 
'oscillons' yj. Due to the subharmonic response of the 
excitation, they resemble during one cycle of the driv- 
ing a circular peak and during the next cycle a circu- 
lar crater. Since the driving is invariant under shifts 
by one period, peaks and craters can exist simultane- 
ously in any given cycle. The excitations are found to 
form bound structures. Thus in the experiments, sta- 
ble pairs, tetramers, and chains of alternating polarity 
are observed. In fact, even the square pattern can be 
thought of as being built up from these individual units 
rather than extended waves. 

Stimulated by these experiments oscillon-type patterns 
have been studied theoretically from a number of differ- 
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ent angles J?1j8p|,[l0| . In fit] oscillons are found within a 
Ginzburg-Landau approach; assuming small amplitudes 
and weak damping of the waves a Ginzburg-Landau equa- 
tion for the wave amplitude is coupled to a mode de- 
scribing the local height of the granular layer. In [|l0| the 
waves are described using a stroboscopic map in which 
the waves arise in a period-doubling bifurcation. A more 
detailed modelling of the dynamics of the granular layer 
is pursued in ||,g,|l^,|ll|,^2| . In |J the granular layer is 
modelled as a sheet that locally follows the dynamics of 
an inelastic ball whereas in j^] the core ingredients are the 
free flight of the granular particles and diffusion. While 
in @,||||[l0) oscillon-type structures have been found nu- 
merically, the mechanism that is responsible for the lo- 
calization of the waves has not been clearly identified. 

In this paper we wish to emphasize two main points. 
First, we demonstrate that oscillon-type patterns are nei- 
ther specific to parametrically excited waves nor to gran- 
ular media. By studying a suitable extension of the Swift- 
Hohenberg equation we show that stable oscillons may 
be expected quite generally in systems that undergo a 
sufficiently hysteretic transition to square patterns. The 
mechanism of localization is provided by non-adiabatic 
effects HIO 

, fnii|]l6f| . In a very recent study of a some- 
what simpler Swift-Hohenberg equation p7[ similar re- 
sults were found. We present here a detailed investiga- 
tion of the range of existence of oscillons and their bound 
states. Second, we compare our general results with the 
experimental observations and comment on what may be 
expected in additional experiments. This will help deter- 
mine whether in this specific system additional mecha- 
nisms related to granular media are important, such as a 
coupling to an additional field J7j. 



II. THE MODEL 

To study the connection between subcritical square 
patterns and oscillons we study a simple Swift- 
Hohenberg-type model. Central for the richness of states 
observed in connection with oscillons is the fact that 
there are two types of oscillons with equal and oppo- 
site polarity. The two types are related by a reflection 
symmetry that arises from the discrete time-translation 
symmetry of the system. Within the framework of a 
Swift-Hohenberg equation this leads to the requirement 
that the equations be equivariant under the reflection 
of the amplitude. We therefore consider the following 
two-dimensional equation for the amplitude ip that char- 
acterizes the pattern, 
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Rip - (<9 2 + lfip + b*p 3 - c4> 5 
+ eV • [(W>) 3 ] - M{V4>) 2 - 



f3 2 tp 2 V 2 ip. (1) 



As in the usual Swift-Hohenberg equation the basic state 
■0 = becomes first unstable to a spatially periodic struc- 
ture with wave number q = 1 at R = 0. The choice of 
nonlinear terms is guided by the following considerations. 
Due to the reflection symmetry ip — > —ip only odd pow- 
ers of tp are possible. In order to obtain a subcritical 
bifurcation the usual Swift-Hohenberg equation has to 
be extended to fifth order. If only the nonlinearities ip 3 
and ip 5 are present no stable square patterns arise. The 
term proportional to e is chosen because it is known to 
favor square patterns (l^Jl^l and the terms proportional 
to the /3i a allow additional tuning possibilities. 

Models of type ([j]) have been used to describe a wide 
range of patterns from Rayleigh-Benard convection to 
wide-area lasers. The connection with the experiments 
on vibrated granular layers can be understood in the fol- 
lowing way. In a weakly nonlinear framework one at- 
tempts to describe the standing waves by an expansion 
based on the mode that destabilizes the basic state with 
undeformed surface, h = 0. The surface deformation of a 
spatially periodic standing wave arising in a subharmonic 
instability can then be written as 



h(r, t) = AH u (t) e Mt /V q ' r + c.c. + h.o.t., 



(2) 



where H u (t) has the same period as the driving, 2tt/lo, 
and A is the small amplitude of the wave. In a Ginzburg- 
Landau approach the amplitude A would be allowed to 
vary slowly in space and time. In order to capture 
the non-adiabatic effects that couple the slow and fast 
space dependence, no slow space dependence can be in- 
troduced. The order parameter tp is therefore allowed 
to vary on a fast spatial scale, which incorporates the 
term e iq r . The fast time dependence, however, is locked 
to the driving and can be kept outside the order pa- 
rameter . From (^) it is apparent that the symmetry 
t — > t + 2ir/uj induces the symmetry ijj(r,t) — » — ip(r, t) of 
the Swift-Hohenberg equation. Note that subharmonic 
instabilities are characterized by a single real Floquct 
multiplier, p, = — 1. Therefore, near onset the order pa- 
rameter -0, describing the excited waves, is real. An ex- 
plicit derivation of a one-dimensional version of (Q) from 
a periodically forced nonlinear Schrodinger equation has 
been performed in the context of phase-sensitive ampli- 
fiers of optical solitons |20j El. In that analysis damping is 
strong and the (strong) driving occurs periodically in a 
shock-like manner. This situation is very similar to that 
encountered in the granular system. 

To obtain an overview of the periodic solutions of ([j]) 
we perform a weakly nonlinear analysis. The ansatz ip = 



'In HJ e = and (d* + l) 2 ip is replaced by {dl - l) 2 ip. 



Ae ix + Be ty + c.c. 
equation 



h.o.t. yields the following amplitude 



d t A = RA - [3 | A | 2 A - 7 | B | 2 A 

-6\ A\ 4 A-2p\ A\ 2 \ B \ 2 A 
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where (3 = -3b + 3e + fa - 3/3 2 , 7 = -66 + 2e + 2/3i - 
6/3 2 , 5 = (-36/3i-9e/3i + 5/3 2 -ll/3i/3 2 + 480cd)/48(i, and 
p = (-3ty?i - ej3i + 3(3f - 7/3i/3 2 + 120cd)/4d. The equa- 
tion for B is obtained by interchanging A and B in (^). 
Thus the bifurcation to stripes and squares is subcritical 
for j3 < and f3 + 7 < 0, respectively. Motivated by the 
experiments, we focus on the regime in which the bifur- 
cation to square patterns is subcritical. Unless otherwise 
noted, the bifurcation to stripes is also subcritical. 



III. OSCILLONS 

To investigate localized, oscillon-type solutions we 
study (Q) numerically using a pseudo-spectral code. A 
single oscillon can be seen in Fig. [j]. Fig. || shows the 
region found numerically for the existence of localized 
states for several values of e for the case 0\ = /? 2 ■ Below 
the solid triangles, oscillons do not exist, whereas above 
the solid circles and squares the oscillons lose stability to 
extended stripes or squares, respectively. Also shown are 
the saddle-node bifurcations for squares and stripes. 
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FIG. 1. Single oscillon for the parameters are 
R = -0.5, b = c = 1.0, e = 0.75, ft = p\ = 2.0. 
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FIG. 2. Numerically obtained saddle-node bifurcations for 
squares (open squares), stripes (open circles), and oscil- 
lons (solid triangle) along with the parameter values at 
which oscillons lose stability to extended squares (solid 
squares) or stripes (solid circles). The parameters are 
b = c= 1.0, fa = jh =2.0. 

We now discuss the mechanism involved in localizing 
these structures. Consider an oscillon as an extremely 
small domain of a patterned state connected by fronts to 
the basic, unpatterned state. Within the framework of 
a Ginzburg-Landau equation, fronts are stationary only 
for a single parameter value. In this setting, the interac- 
tion of fronts is attractive and such a localized state is 
unstable (e.g., pl|,p2|). In |13| the pinning of interfacial 



fronts by the interaction with the underlying small-scale 
periodic pattern has been pointed out. The pinning ef- 
fect is examined in a more quantitative analysis in p4[ . 
It is found that the front is locked for a range of the con- 
trol parameter which is found to scale as e~ a /^'. This 
same pinning was used to explain the existence of local- 
ized pulses in a 1-D quintic Swift-Hohenberg equation 
fusf . We expect the same mechanism to operate in two 
spatial dimensions. Hence, we also expect the range of 
existence of the oscillons to scale in a similar manner. 

To confirm that the pinning is responsible for the lo- 
calization, we investigate the dependence of the existence 
range of oscillons as a function of c. Increasing the prefac- 
tor, c, of the quintic term causes both squares and stripes 
to become less subcritical. Fig. ||a shows a phase diagram 
for the range of existence of the localized structures along 
with the numerical and analytical saddle-node curves for 
squares. Since squares become less subcritical, an in- 
crease in c necessitates an increase in R for finding oscil- 
lons. The oscillons have oscillatory tails that decay expo- 
nentially outward. As R is increased, the spatial decay of 
these tails becomes more gradual producing larger radial 
oscillations. Hence, the front locking becomes weaker and 
the localized solution loses stability to a periodic pattern. 
To demonstrate the scaling of the locking range, we plot 
the range in the control parameter R where a single os- 
cillon exists as a function of the inverse square root of 
the midpoint of this range. The log-linear plot in Fig. Bp 
shows that the range does scale proportional to e" a ^ R . 
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FIG. 3. (a) Oscillons exist between the solid triangles and 
the open triangles. The numerical saddle-node bifurcation 
for squares (open squares) along with the analytical sad- 
dle-node curve (thin solid line) are shown. The thick solid 
line shows the analytical curve where AE = (cf. §)). (b) 
The existence range of a single oscillon is plotted as a func- 
tion of the inverse square root of t he m idpoint of the range. 
The dashed line is a fit to ne~ a ^ Rm . The parameters are 
b= 1.0, e = 0.75, /3i = (3 2 = 2.0. 

For (3i — P2 — 2/ equation (|l]) possesses a Lyapunov 
Functional £, 



bip 4 + -cip 6 
6 



C = J &{-ii^ 2 + i[(v 2 + l)Vf - \ 

- ie(V^) 4 + i/^V 2 (^ 2 )} ee j dxE(x) 



(4) 



with ip t = —6£/5ijj. In this case the 'energy' difference 
AE between different states of the system governs the 
motion of a single front. For the stable zero state E is 
0, while for squares it is -RA 2 /2 + (j3 + "f)A A /A + (6 + 
3p)A 6 /6 within the weakly nonlinear theory (|J). The 
values of R at which both states have the same energy 
(AE = 0) and at which a single front would be station- 
ary is shown as the thick curve in Fig. ||a. The existence 
range of oscillons lies to the right of this curve, corre- 
sponding to parameters for which the square pattern has 
a lower energy. We attribute this difference to the at- 
tractive interaction between the fronts making up the 
oscillon. 

When the oscillon loses stability to an extended pat- 
tern state, stripes or squares emerge depending on other 
parameters. Fig. || shows intermediate stages for a sin- 
gle oscillon in the center evolving to stripes and squares, 
respectively. As can be seen in the transition to stripes, 
the oscillon tends to elongate in one direction and gener- 
ate small stripes in the other direction, eventually filling 
the entire space with stripes. Note that the square pat- 
tern, although being unstable to stripes, still exists in 
this regime (cf. Fig. |^). The transition to squares is gen- 
erally characterized by more oscillons appearing on the 
perimeter of the initial oscillon. These features seen in 
the numerical simulations are similar to the transitions 
observed in the experimental granular system. 
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FIG. 4. The transition to stripes (a) and squares (b) start- 
ing from a single oscillon in the center. The parameters are 
b = c = 1.0, pi = = 2.0 and for (a) R = -0.264, e = 0.5 
and (b) R = -0.205, e = 0.75. 

To investigate whether the general locking mechanism 
is sufficient to explain the experiments or whether as- 
pects have to be taken into account that are specific to 
the granular system, we compare our findings with the 
experimental results in more detail. The experimental 
phase diagram in Fig. || jj| shows that oscillons exist only 
below the patterns for a mid-range of frequencies. In- 
creasing the driving amplitude T, patterns appear at the 
solid diamonds. For decreasing T, extended patterns dis- 
appear at the solid circles to become either the flat layer 
or oscillons which exist all the way to the solid squares. 
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FIG. 5. The experimental phase diagram JO]. For increas- 
ing r, patterns appear at the solid diamonds. As F is de- 
creased the extended patterns exist until the solid circles be- 
yond which oscillons exist in the mid-frequency range until 
the solid squares. 

Within the minimal framework presented here, oscil- 
lons should exist where squares are strongly subcritical. 
However, in the low frequency range where squares ex- 
hibit strong hysteresis, oscillons are not seen experimen- 
tally. We propose the following explanation. In the ex- 
periments, the driving amplitude is increased into the 
pattern regime which results in a pattern with wavenum- 
ber q « q c . Then the amplitude is decreased and oscillons 
emerge depending on where the pattern with this wave 
number goes unstable. Fig. || shows the neutral stabil- 
ity curve for the basic state -0 = and the numerically 
determined range over which square patterns (restricted 
to one wavelength) are stable for two sets of parameter 



values. Decreasing R below the value marked by solid 
squares, the square pattern goes to zero in a saddle node 
bifurcation. The open squares indicate the value of R 
where the square pattern seems to undergo an instability 
involving the wavelength of the pattern. Note that we 
do not associate a wavenumber q with the oscillon and 
the shaded region refers only to the R- values. As the 
control parameter R is decreased, squares remain in ex- 
istence until the solid saddle-node points, although they 
may lose stability earlier. In Fig. ga the oscillon range 
extends below the saddle-node for squares, so that oscil- 
lons are expected as R is decreased below the existence of 
squares. In contrast, Fig. ^|b shows a situation in which 
the oscillon range is above these points near q c = 1, and 
we do not expect to obtain oscillons by decreasing R fur- 
ther below the saddle-node. 
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FIG. 6. The neutral stability curve for the basic state and 
the numerical saddle-nodes for squares (solid squares) are 
shown as a functions of the wavenumber q. Open squares 
denote cases where the single square loses stability. For 
the parameters b = c = 1.0, e = 4.0, /?i = 0.0, /?2 = 2.3 
(a), the oscillons exist for values of the control parame- 
ter below the saddle-node for squares. For the parameters 
b = c = 1.0, e = .75, /3i = /3a = 2.0 (b), the oscillons no 
longer exist before the saddle-node for squares is reached by 
decreasing R. 

Fig. ^ shows the bifurcation diagram obtained from 
numerical simulations of (hi) for the same parameters as 
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in Fig. |6|a. We note that stripes are supercritical in this 
case. Here we see that a single oscillon exists for a range 
of R values below the saddle-node for the square pattern 
with wavenumber q = 1. 

1.75 | , 1 , 1 , 1 



be 'embedded' in the range of square existence as sug- 
gested by Fig. IJd. In that case the oscillons would not 
be accessible by changing the amplitude of forcing. We 
suggest that oscillons may be accessible in this regime by 
starting with an oscillon in the mid-frequency range and 
then decreasing the frequency. 
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FIG. 7. Bifurcation curves obtained for a single oscil- 
lon (solid triangles) and squares with q = 1.0 (solid 
squares) and q = 0.875 (open squares). Parameters are 
b = c = 1.0, e = 4.0, /3i = 0.0, & = 2.3 (cf. Fig. |a). 

Starting with a slightly perturbed, extended square 
pattern and then decreasing the control parameter below 
the saddle-node, we find that the pattern breaks down 
and some oscillons persist. Figs. |a-c show a case in 
which several oscillons survive, while Figs. ||d-f show a 
case in which only one oscillon remains. These results 
are consistent with the observation of oscillons existing 
below the periodic patterns in the granular experiments. 
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FIG. 8. (a) An initially perturbed square pattern with 
R — —0.45, then R is decreased to —0.48 (b) where parts 
of the pattern decay (t = 15), and (c) several oscillons persist 
(t — 1000). (d) Again, the perturbed square pattern with 
R = —0.45, (e) followed by rapid decay of the pattern for 
R — —0.504 (t = 10), (c) where only one oscillon survives 
(t — 2000). All other parameters as in Fig 

In light of the experiments for lower frequencies, the 
oscillons may still exist, but their range of existence may 




IV. INTERACTION OF OSCILLONS 

Since the oscillons have oscillatory tails that decay ex- 
ponentially outward, one expects them to interact and 
form bound states. Fig. || shows several of these bound 
states found in the numerical simulations which are simi- 
lar to states observed in the experimental system. These 
include dipoles, tetramers, and chains of alternating po- 
larity. 




FIG. 9. Bound states of opposite polarity: (a) dipole, (b) 
tetramer, and (c) chain of 5 alternating oscillons for the pa- 
rameters as in Fig. [j]. 

In addition to seeing oscillons of opposite polarity bind- 
ing, we also find bound states of like polarity, which have 
not been seen in the experiments. The separation dis- 
tance between like pairs is greater than for opposite pairs. 
As can be seen in the cross-section of a like pair, Fig. [l(], 
the spatial oscillations between the pair lock the two os- 
cillons together, but do not create a third negative oscil- 
lon between them as can be seen by comparison with the 
H 1- chain (dashed line). However, increasing the pa- 
rameter R causes the oscillations to decay slower in space. 
Eventually a value of R is reached where the minimum 
exceeds a threshold beyond which the amplitude grows 
to form an oscillon of opposity polarity and becomes a 
H h chain. 
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FIG. 10. (a) Two bound oscillons of like polarity and (b) 
their cross-section (solid line) along with the cross-section for 
a chain of alternating polarity (dashed line). For the pair 
R = -0.470 while for the chain R = -0.469. All other pa- 
rameters as in Fig. jjj. 

Fig. [ll] shows the existence range of ++ structures for 
two values of c (dark region). The border between the 

dark and grey regions marks the transition to a H h 

structure. To the right of the grey region the H h chain 

evolves to a periodic pattern. This result shows that the 
ratio of the like-pair range to the single-oscillon range is 
greatly reduced by slightly increasing c. Increasing c fur- 
ther makes it difficult to find a stable ++ structure. We 
expect that a bound pair of like polarity may be simi- 
larly difficult to find in experiments because they exist 
only over a small range. Conversely their range of exis- 
tence may be larger in regions of stronger subcriticality. 
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FIG. 11. Existence range for a pair of like oscillons is shown 
in black. The left side of the grey region marks the transition 
to a chain of three oscillons of alternating polarity (cf. Fig. 

|l(i|b). The H \- chain persists up to the right end of the grey 

area. Parameters are as in Fig. 

Since the peak-peak distance for two like pulses is 
greater than that for two pulses of opposing amplitude, 
we expect the binding between the ++ state to be much 
weaker than the binding of the — H state. In order to 
study the 'strength' of the binding, we add bounded noise 



to the solution at each time step and measure the peak- 
peak distance. The results for several noise strengths 
are shown in Fig. 03. Fig. 12 1 shows that the opposing 



pair remains bound for all three noise levels. In contrast, 
Fig. |l2| b shows that the like pair is more susceptible to 
the same values of noise. In fact, for the mid-level noise, 
the like pair creates a third negative oscillon in the cen- 
ter and for the highest-level noise, the pair quickly loses 
stability and by 40000 time steps has become two perpen- 
dicular chains of alternating polarity. The effect of noise 
in the experimental system combined with the small ex- 
istence range for oscillons of like polarity may contribute 
to them not being observed in the experiments. 



4.5 



a 



■a 

a 
p 

- 

a 
& 

0J 

[B 

*6 

Ct 

a 



4.0 



3.5 



10000 20000 30000 40000 
time steps (dt = 0.005) 



7.5 



u 
s 

es 



S 

o 



transition 
to+-+ 



•3 7.0 



C3 
■- 

es 
a. 

in 
C3 

aj 
& 




++ 



6.5 



10000 20000 30000 40000 
time steps (dt = 0.005) 

FIG. 12. (a) The peak-peak distance for a pair of oscil- 
lons of opposite polarity is plotted for the unperturbed pair 
(straight line) and for three added noise levels, 0.075 (bot- 
tom), 0.15 (middle), 0.3 (top, grey), (b) The same is plotted 
for a like pair. The parameters are as in Fig. ^. 

The interaction of like pairs can also result in forming 
triad structures of two like oscillons bound to a third op- 
posite oscillon in a triangular pattern as shown in Fig. |l3|. 
The triad breaks the spatial reflection symmetry and in 
this non-potential case the triad travels in the direction 
of the two leading oscillons. 
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FIG. 13. Triad propagating downward for 
R = -1.5, b = c = 1.0, e = (3 2 = 4.0, A = 0.0. 

V. CONCLUSION 

We have presented the results of investigating a Swift- 
Hohenberg model for the existence of two-dimensional 
localized oscillons. Oscillon-type structures are a gen- 
eral feature to be expected in systems where the tran- 
sition to squares is sufficiently subcritical so as to allow 
localization caused by the interaction of the large-scale 
envelope with the small-scale underlying periodic pat- 
tern. Hence the observance of oscillons should not be 
restricted to granular materials. For example, we would 
also expect them to be observed in Faraday experiments 
of shear-thinning fluids for which the viscosity decreases 
with the flow, which should make the transition subcriti- 
cal. Since the subharmonic response results in the reflec- 
tion symmetry tp — > — -0, one obtains equal oscillons of 
opposite polarity. For transcritical hexagons similar lo- 
calized structures are also obtained, but since there is no 
reflection symmetry, one does not obtain two different, 
equivalent oscillon-structures pj| . 

A comparison with the granular experiments will show 
whether this general mechanism is sufficient to capture 
the observed oscillon. Experimental tests to determine 
the validity of this minimal framework involve the range 
of existence of oscillons and of like-pairs in the low fre- 
quency range. We expect that oscillons could be observed 
experimentally in the low frequency regime by chang- 
ing the experimental protocol. We also note that bound 
states of like-polarity are shown numerically to be not as 
robust as those of opposite polarity. Their range of ex- 
istence may be larger in regions of greater subcriticality, 
e.g. lower frequencies, and their observance more likely 
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